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SUMMARY 


This  Final  Report  describes  the  research  work  on  the  contract  of  investigation  entitled  by 
Study  of  Cryogenic  Plasma  in  Superfluid  Liquid  Helium  (Principal  investigator:  O.  Ishihara, 
AOARD-04-4018)  for  the  period  from  June  7,  2004  to  August  6,  2005.  The  research  has  been 
carried  out  at  Yokohama  National  University  in  collaboration  with  Prof.  K.  Minami  of  Tokyo 
Denki  University,  a  research  associate  Dr.  Masako  Shindo,  a  doctoral  student  Mr.  Chikara  Kijima, 
and  two  master  students.  A  preliminary  experiment  to  produce  a  DC  discharge  plasma  in  a  vapor 
of  liquid  helium  started  and  a  complex  plasma  with  charged  fine  particles  in  the  cryogenic 
environment  has  been  produced  successfully.  The  analytical  study  on  the  experimental  results 
completed  at  Niigata  University  for  the  pulsed  discharge  in  a  microwave  cavity  in  liquid  helium 
revealed  the  nature  of  decaying  plasma  in  liquid  helium. 


I.  Objectives 

Our  overall  goal  is  to  study  the  cryogenic  plasma  experimentally  and  theoretically  and 
to  reveal  novel  natures  of  cryogenic  plasma  produced  by  a  pulsed  or  a  DC  discharge  above  or  in 
superfluid  liquid  helium . 

II.  Status  of  the  Research 

(1)  Preparation  of  the  Dewar  Bottles 

The  Physics  Department  in  the  faculty  of  engineering  at  Yokohama  National 
University  was  established  in  1998  and  the  principal  investigator  Osamu  Ishihara 
joined  the  faculty  of  Yokohama  National  University  in  June  of  1999.  The  plasma  lab 
was  established  in  the  August  of  2000  immediately  after  the  completion  of  the  new 
departmental  building.  Initially,  major  efforts  in  the  research  were  devoted  to  the 
theoretical  and  computational  plasma  physics.  The  collaboration  with  Prof.  K. 
Minami  of  Niigata  University  started  in  2000  and  the  AOARD  funding  started  in 
2002  (AOARD -00- 40 14,  Study  of  cryogenic  plasma  produced  by  a  pulsed  discharge  in 
Superfluid  liquid  helium).  The  plasma  lab  expanded  in  space  and  the  complex 
plasma  experimental  lab  in  the  Department  of  Physics  at  Yokohama  National 
University  was  established  in  March  2004. 

The  principal  investigator,  Professor  Osamu  Ishihara,  has  been  conducting  the  research  in 
collaboration  with  Dr.  Masako  Shindo,  Research  Associate  joined  in  the  group  in  April  of  2002, 
and  a  doctoral  student  Mr.  Chikara  Kojima  who  joined  in  the  group  in  April  of  2004  and  two 
master  students,  Mr.  M.  Kugue  and  Mr.  T.  Maesawa.  Professor  Kazuo  Minami,  retired  from 
Niigata  University  and  now  with  Tokyo  Denki  University,  is  consulting  us  on  the  research. 

Following  the  retirement  of  Professor  Minami  from  Niigata  University  in  March,  2004, 
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a  cryostat  system  including  a  Dewar  vessel  was  transferred  to  Yokohama  National  University.  A 
preliminary  DC  discharge  in  the  vapor  of  liquid  helium  in  the  Dewar  bottle  was  performed  and 
preliminary  probe  measurements  have  been  carried  out.  The  setup  is  schematically  shown  in 
Fig.  1,  where  two  thin- walled  bottles  nested  one  inside  the  other  and  the  outer  bottle  is  sealed  at 
the  ends.  The  one  end  of  the  inner  bottle  needs  to  be  pumped  to  maintain  the  vacuum  for 
insulation.  The  DC  discharge  in  a  vapor  of  liquid  helium  was  accomplished  by  applying  about 
400V  between  a  meshed  stainless  anode  and  a  stainless  disk  cathode  separated  by  6  mm.  The 
plasma  was  characterized  by  the  electron  temperature  about  0.3  eV  and  the  plasma  density  about 
4x  107  cm"3.  We  note  that  the  temperature  of  He  vapor  in  the  Dewar  increased  up  to  86K  during 
the  discharge. 
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Fig.l  Schematic  diagram  of  the  experimental  apparatus.  A  Dewar  bottle  has 
two  thin-walled  bottles  nested  one  inside  the  other  and  the  outer  bottle  is 
sealed  at  the  ends.  The  electrodes  are  set  above  the  liquid  helium  (LHe). 
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A  larger  Dewar  vessel  to  accommodate  a  DC  discharge  chamber  was  designed,  ordered 
and  delivered.  It  is  aimed  to  produce  a  plasma  in  a  vapor  of  liquid  helium  in  a  confined  chamber 
which  is  immersed  in  the  liquid  helium  rather  than  in  a  vapor  of  liquid  helium.  Since  the 
delivered  vessel  had  a  defect,  the  vessel  had  to  be  repaired  at  the  JEC-Tori  Corporation  of 
Saitama  Prefecture.  While  the  Dewar  vessel  used  in  Niigata  is  rather  small,  10  cm  in  inner 
diameter,  the  new  Dewar  vessel  is  16cm  in  inner  diameter  to  accommodate  a  DC  discharge 
chamber  of  10  cm  in  diameter.  The  Dewar  vessel,  as  shown  in  Fig.  2,  has  two  thin-walled  bottles 
nested  one  inside  the  other  and  both  bottles  are  sealed  at  the  ends.  A  DC  discharge  chamber,  to  be 
immersed  in  the  liquid  helium,  was  designed  ordered  and  delivered,  while  the  original  chamber 
delivered  needed  to  get  fixed  due  to  a  defect  at  an  injection  valve  found  later.  The  JEC  Tori 
Corporation  repaired  the  discharge  chamber  and  delivered  the  chamber  to  us. 
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Fig.  2  A  larger  Dewar  vessel,  designed  and  constructed  at  Yokohama 
National  University  with  a  DC  discharge  chamber. 
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(2)  Discharge  in  a  chamber  placed  in  the  Dewar  vessel 

A  preliminary  discharge  between  a  brass  rod  cathode  and  a  brass  rod  anode  in  the  DC 
discharge  chamber,  placed  in  the  Dewar  vessel,  produced  a  plasma  in  a  helium  gas  of  0.5  Torr 
maintained  to  be  77K  by  the  heat  exchange  through  the  chamber  wall.  The  electrodes,  both 
10mm  in  diameter  and  40mm  in  length,  were  placed  10  mm  apart  in  the  chamber  and  were 
applied  with  the  DC  voltage  of  about  600  V.  When  the  chamber  was  immersed  in  liquid  nitrogen 
in  the  Dewar  vessel,  liquid  nitrogen  instead  of  liquid  helium  at  a  preliminary  stage,  the  produced 
plasma  was  rather  unstable  to  maintain  the  constant  density.  A  Langmuir  probe  with  0.5  mm  in 
diameter  and  10  mm  in  length  measured  the  electron  temperature  of  order  1  eV  and  revealed  the 
slight  decrease  in  electron  temperature  in  the  presence  of  liquid  nitrogen  around  the  DC  chamber 
compared  with  the  electron  temperature  without  liquid  nitrogen.  While  we  introduced  liquid 
helium  into  the  DC  chamber,  a  leakage  was  found  at  a  valve  attached  at  the  top  of  the  DC 
chamber  to  introduce  liquid  helium.  Our  efforts  to  produce  a  stable  plasma  in  a  vapor  of  liquid 
helium  continue  and  the  plasma  parameters  have  been  measured  yet  with  insufficient  accuracy. 
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Fig.  3  DC  discharge  chamber  to  accommodate  liquid  helium  (LHe)  and  to  be  immersed 
in  liquid  helium  in  the  Dewar. 
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(3)  Production  of  a  steady  state  plasma  in  a  vapor  of  the  liquid  helium 

A  cryogenic  plasma  in  a  steady  state  has  been  established  by  an  electrical  discharge 
(AClOkHz,  ~lkV)  in  a  vapor  (2  Torr)  above  the  liquid  helium  surface.  The  electrical 
discharge  in  the  vapor  necessarily  brings  a  heat  in  the  environment  resulting  in  the  boiling  of 
the  liquid  helium.  To  overcome  the  rise  of  the  surrounding  plasma  temperature  as  well  as  the 
instability  of  the  surface  of  the  liquid  helium,  a  special  purpose  discharge  unit  has  been 
developed.  The  unit  is  characterized  by  a  region  of  plasma  production  isolated  from  the 
surroundings  by  a  small  glass  Dewar.  Two  identical  plane  discharge  electrodes  have  coaxial 
holes  at  the  center  through  which  the  helium  gas  was  fed  from  the  upper  hole.  A  plasma  is 
produced  by  a  discharge  between  the  electrodes  and  the  gas  flow  takes  the  plasma  outside  of 
the  discharge  region  through  the  bottom  hole.  The  rise  of  the  plasma  temperature  once  outside 
of  the  small  glass  Dewar  is  successfully  suppressed  and  the  plasma  temperature  decreases 
toward  the  surface  of  the  liquid  helium.  The  cold  plasma  in  a  steady  state  is  thus  successfully 
produced.  The  first  report  on  this  subject  will  be  presented  by  a  doctoral  student  Mr.  Chikara 
Kojima  at  the  domestic  meeting  in  September  (see  Appendix  3). 

(4)  Complex  plasma  in  the  ultra-cold  environment 

A  complex  plasma  is  characterized  by  a  group  of  negatively  charged  macro  particles 
embedded  in  a  plasma.  The  group  of  fine  particles  shows  collective  behaviors,  while  plasma 
particles  show  collective  behavior  independently  with  a  time  scale  much  faster  and  a  space 
scale  much  shorter  than  those  for  the  group  of  fine  particles.  The  complex  system  including 
fine  particles  in  a  plasma  is  characterized  by  the  unique  interaction  of  the  two  collective 
modes.  The  fine  particles  (often  called  dust  particles)  are  strongly  coupled  through  Coulomb 
interaction  with  a  large  coupling  constant  T  defined  by  the  ratio  of  the  Coulomb  energy  and 
the  thermal  energy  of  particles: 

^  _  Coulomb  Energy  _  (Zdef 
Thermal  Energy  4 n£0ricTd  ’ 

where  Zde  is  the  charge  on  the  surface  of  a  fine  particle,  r  is  the  average  inter-particle 
distance,  k  is  the  Boltzmann  constant,  and  Td  is  the  temperature  of  fine  particles.  We 
used  here  the  subscript  d  for  dust  particles  (fine  particles).  Ordinary  plasmas  with  electrons 
and  ions  are  characterized  by  rather  a  small  coupling  constant,  T « 1 ,  since  the  kinetic 
energy  of  plasma  particles  are  dominant  over  the  Coulomb  energy  allowing  the  random 
motion  of  plasma  particles.  The  coupling  constant  for  fine  particles  could  be  much  larger  than 
1 ,  either  through  \Zd\ » 1  or  KTd  —>  0 . 

In  an  extremely  low  temperature  environment,  the  thermal  energy  of  the  particles 
are  suppressed  resulting  in  a  large  coupling  constant.  We  have  injected  fine  acrylic  particles 
of  5pm  in  diameter  in  a  vapor  above  liquid  helium  at  about  IK.  The  pressure  measured  at  the 
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upper  part  of  the  Dewar  was  about  2  Torr.  The  discharge  electrodes  were  set  200mm  above 
the  surface  of  liquid  helium,  and  the  distance  between  electrodes  was  5mm.  The  10kHz  AC 
voltage  of  lkV  was  applied  to  the  driving  electrodes.  The  surrounding  of  electrodes  were 
covered  by  the  insulator  to  suppress  the  rapid  rise  of  the  temperature  of  both  liquid  helium  and 
the  vapor  gas.  The  produced  plasma  was  diffused  toward  Liquid  helium  surface  through  a 
5mm  hole  in  lower  electrode.  It  was  confirmed  that  an  electron  temperature  in  the  diffusion 
region  decreased  due  to  the  presence  of  liquid  helium  by  double  probe  measurements.  The 
acrylic  particles  were  dropped  from  the  upper  part  into  plasma  through  the  5mm  holes  at  the 
center  of  the  electrodes.  They  are  expected  to  be  charged  negatively  and  blew  out  into  the 
diffusion  region.  The  diffusion  region  was  covered  with  an  acrylic  tube  of  50mm  in  diameter 
to  confine  the  particle  in  the  radial  direction,  and  a  negatively-biased  plate  electrode  of  50mm 
in  diameter  was  set  100mm  under  the  discharge  electrodes.  The  He-Ne  laser  light  was 
irradiated  from  outside  though  a  10mm- width  slit  on  the  glass  Dewar,  which  was  silver  plated 
except  for  the  slit  area,  and  the  acrylic  particles  float  10mm  above  the  plate  electrodes  in  a 
plasma.  The  preliminary  results  will  be  presented  in  the  coming  Fall  meeting  of  the  Physical 
Society  of  Japan  by  a  research  associate  Dr.  Masako  Shindo  (See  Appendix  3). 

(5)  Study  of  the  equation  of  state 

We  have  shown  that  the  long-range  interaction  energy  between  fine  particles  placed 
in  the  plasma  exhibits  a  unique  feature  with  a  repulsive  nature  for  short  distances  and  with  an 
attractive  nature  for  larger  distances  as  shown  in  Fig.  4  (O.  Ishihara  and  N.  Sato,  Attractive 
Force  on  Like  Charges  in  a  Complex  Plasma,  Physics  of  Plasmas,  12,  070705  (2005),  see 
Appendix  4).  Such  a  finding  suggests  the  possibility  of  the  phase  transition  in  the  complex 
plasma.  Although  the  topic  of  the  equation  of  the  state  was  not  originally  included  in  the 
proposal,  the  study  of  dusty  plasmas  prompted  us  to  direct  our  attention  to  the  equation  of 
state  applied  to  a  complex  plasma  in  a  cryogenic  environment.  The  long  range  interaction 
energy  was  formulated  by  the  consideration  of  free  energy  available  in  the  complex  plasma 
where  the  charge  neutrality  is  satisfied  by  the  presence  of  electrons,  ions,  and  the  charged  fine 
particles.  With  the  explicit  formulation  of  the  interaction  energy,  the  second  Virial  coefficient 
will  be  calculated.  And  the  equation  of  state  will  give  us  a  clue  to  the  study  of  the  critical 
phenomenon  in  the  complex  plasma.  The  possible  presence  of  a  critical  point  in  the  complex 
plasma  in  the  cryogenic  environment  is  still  an  open  question  both  theoretically  and 
experimentally. 


Fig.  4  Interaction  energy  of  fine 
particles  in  a  complex  plasma. 


(6)  Analysis  of  the  ultracold  plasma  produced  by  a  pulse  discharge 

After  all  the  experiments  at  Niigata  University  completed,  the  data  taken  for  the 
pulsed  discharge  were  analyzed  through  the  initiative  of  Professor  Minami.  Our  focus  was  on 
the  understanding  of  the  decaying  plasma  produced  by  a  high-voltage  (~20kV)  pulse  of 
duration  on  the  order  of  micro  seconds  in  an  X-band  mode  cylindrical  cavity  filled  with 
Liquid  helium.  Although  the  transmission  signals  through  the  cavity  show  clearly  the  presence 
of  decaying  plasma  in  liquid  helium  with  decay  time  on  the  order  of  a  few  hundred  ms  for 
temperature  of  a  few  K  at  saturating  vapor  pressure,  our  estimated  mass  of  charged  carriers 
seems  to  be  much  less  than  expected  from  the  earlier  theoretical  predictions.  The  paper  is  now 
in  print  (see  Appendix  1)  and  the  results  were  presented  at  the  international  conference  (see 
Appendix  2). 


III.  Publications/  Presentations 

1.  Publication  (Appendix  1) 

K.  Minami,  C.  Kojima,  and  O.  Ishihara,  Microwave  measurement  of  Decaying  Plasma  in 
Liquid  Helium,  IEEE  Transaction  on  Plasma  Science,  August  issue  (2005). 

2.  International  Conference  Proceeding  (Appendix  2) 

C.  Kojima,  O.  Ishihara,  and  K.  Minami,  Study  of  Decaying  Plasma  in  Liquid  Helium,  19th 
International  Conference  on  Numerical  Simulation  of  Plasmas  and  7th  Asian  Pacific 
Plasma  Theory  Conference  (July  12-15,  2005,  Nara,  Japan)  PI -65,  pp.  15 1-152. 

3.  Domestic  Conference  Abstracts  (Appendix  3) 

M.  Shindo,  C.  Kojima,  O.  Ishihara,  Observation  of  dust  particles  in  a  plasma  under  cryogenic 
environment,  Fall  Meeting,  Physical  Society  of  Japan  (2005.9.19-22,  Kyoto)  21  pWG- 
1 . 

C.  Kojima,  M.  Kugue,  T.  Maezawa,  M.  Shindo  and  O.  Ishihara,  A  plasma  production  in  a 
region  of  low  temperature,  Fall  Meeting,  Physical  Society  of  Japan  (2005.9.19-22,  Kyoto) 
22  aWG-2. 
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IV.  AOARD  site  visits 

August  6,  2004 


AOARD  program  manager  Dr.  Misoon  Y.  Mah 
September  24,  2004 

AOARD  program  manager  Dr.  Misoon  Y.  Mah, 

AOARD  Science  Advisor  Dr.  Takeo  Miyazaki, 

Air  Force  Research  Laboratory  (Kirtland  AFB)  Senior  Research  Physicist  Dr.  S.  Joe  Yakura 
AF  Research  Laboratory  (Kirtland  AFB)  Chief  of  Applications  Branch  Maj.  Jeffrey  C.  Wiemeri 


Maj.  Wiemeri,  Dr.  Yakura,  Dr.  Mah,  Dr.  Miyazaki,  Prof.  Ishihara 
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V.  Lab  photos 


Osamu  Ishihara,  Masako  Shindo  and  three  graduate  students 


Professor  Ishihara’s  lab  was  introduced  in  an  annual  campus  guide  YNU  2005. 
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APPENDIX  1 
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IEEE  TRANSACTIONS  ON  PLASMA  SCIENCE  1 

Microwave  Measurement  of  Decaying 
Plasma  in  Liquid  Helium 

Kazuo  Minami,  Member,  IEEE #  Chikara  Kojima,  Takeo  Ohira,  and  Osamu  Ishihara,  Fellow ;  IEEE 


Abstract — Decaying  plasma  in  a  microwave  cavity  filled  with 
liquid  helium  (LHe,  hereafter)  is  studied.  An  X-band  TEn3  mode 
cylindrical  cavity  is  filled  with  LHe,  and  a  high-voltage  pulse  with 
duration  7  /is,  voltage  V  <  20  kV,  and  current  I  <  360  A 
is  applied  between  a  tungsten  needle  electrode  and  a  thin  stain¬ 
less-steel  mesh  which  separates  an  adjacent  small  discharge  space 
from  the  cavity.  The  transmission  signals  through  the  cavity  show 
the  presence  of  decaying  plasma  in  LHe  with  decay  time  on  the 
order  of  200  ms  for  a  temperature  range  from  4.2  to  23  K  at 
saturating  vapor  pressure.  The  response  signals  are  suppressed  in 
the  pressurized  LHe,  as  well  as  in  superfluid  LHe  below  2.17  K. 
Our  estimated  mass  of  the  charged  carriers  just  after  the  pulsed 
discharge,  assuming  a  diffusion  loss  mechanism,  is  much  less 
than  the  value  known  as  an  effective  mass  of  heavy  ions  which 
were  produced  by  a  steady-state  ion  source  in  LHe  measured  by 
Poitrenaud  and  Williams  (1972).  Our  experimental  observation 
suggests  that  the  heavy  effective  mass  of  ions  in  LHe  might  take 
time  to  develop  to  its  hill  extent. 

Index  Terms — Cavity,  cryogenic  plasma,  diffusion,  electron  bub¬ 
bles,  ion  clusters,  liquid  helium,  microwave,  pulsed  discharge,  re¬ 
combination,  superfluid. 


I.  Introduction 

ULTRACOLD  plasma  is  a  new  subject  of  research  in  plasma 
science  that  is  currently  in  rapid  progress.  While  a  new 
technology  including  laser  cooling  has  been  utilized  to  pro¬ 
duce  an  ultracold  plasma  [11,  we  have  successfully  produced 
cryogenic  decaying  plasma  in  liquid  helium  (LHe,  hereafter), 
as  well  as  in  cryogenic  gas  helium  with  temperature  below  4.2 
K  by  means  of  a  pulsed  discharge  [2]-[5].  Cryogenic  plasma 
in  LHe  is  known  to  be  characterized  by  features  such  as  elec¬ 
tron  bubbles  and  positive  ion  clusters  [6]— [8],  In  our  earlier 
experiments,  we  succeeded  in  producing  a  high-density  local¬ 
ized  plasma  in  LHe.  High-voltage  pulsed  discharge  with  needle 
electrodes  in  LHe  produced  the  localized  plasma  with  density 
more  than  101*  r*m“3  and  electron  and  ion  temperatures  of  a 
few  electronvolts  [2],  [3].  Upon  termination  of  the  discharge, 
the  plasma  remained  as  long  as  on  the  order  of  10  /is  in  LHe. 
Some  of  ions  and  electrons  remained  longer  in  LHe  before  they 
recombine  [5].  Such  a  cryogenic  plasma,  or  positive  and  nega¬ 
tive  ion  (p-ion  and  n-ion,  hereafter)  mixture  dissolved  in  LHe, 
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is  the  subject  of  our  interest  in  this  paper.  In  our  previous  paper 
[4],  a  plasma  in  cryogenic  helium  gas  in  a  TEnu  mode  S-band 
cavity  with  very  long  decay  time  on  the  order  of  1  s  was  ob¬ 
served.  In  this  paper,  we  apply  the  microwave  cavity  method  to 
measure  decaying  plasma  in  LHe.  Late  afterglow  plasma  in  an 
X-band  cavity  filled  with  LHe  in  the  temperature  range  from  4.2 
to  1.69  K  is  studied.  Slow  decaying  plasmas  in  LHe  detected  as 
long  as  200  ms  after  the  discharge  are  reported  in  this  paper.  The 
plasma  is  observed  in  LHe  at  a  saturating  vapor  pressure  with 
temperature  above  2.3  K.  Although  the  cavity  method  may  not 
be  fully  applicable  to  the  present  study  because  of  the  nature  of 
the  local  plasma  production,  our  estimated  density  based  on  the 
cavity  method  suggests  an  effective  mass  of  ions  much  different 
from  that  expected  by  the  conventional  theory.  A  simple  diffu¬ 
sion  loss  model  for  a  temporal  evolution  of  the  localized  plasma 
produced  by  the  pulsed  discharge  in  LHe  reveals  the  presence  of 
charged  carriers  with  masses  much  lighter  than  those  expected 
in  LHe. 

The  organization  of  the  paper  is  as  follows.  In  Section  II,  the 
principle  of  cavity  measurements  that  were  used  in  our  previous 
paper  [4]  are  reviewed  briefly.  The  experimental  setup,  proce¬ 
dure,  and  experimental  results  are  described  in  Section  III.  in 
Section  IV,  expected  temporal  evolution  of  localized  decaying 
plasma  in  LHe  is  analyzed  to  compare  with  the  experimental  re¬ 
sults.  Discussion  and  conclusion  are  given  in  Section  V. 


II.  Principle  of  Cavity  Measurements 

In  our  previous  paper  on  cryogenic  gas  plasma  [4],  cavity 
measurement  was  applied  to  detect  the  density  n  in  the  de¬ 
caying  plasma.  The  method  is  now  extended  to  the  measure¬ 
ment  of  localized  plasma  in  LHe  in  an  X-band  cavity.  The  an¬ 
gular  frequency  u,*  of  the  sensing  microwave  is  adjusted  to  be 
-e  =  u,*o  +  Bo,  as  shown  in  Fig.  1  in  [4],  where  u/o  and  B0  arc 
the  resonant  angular  frequency  of  the  cavity  without  plasma  and 
half  half- width  of  the  resonant  curve,  as  shown  in  Fig.  1  in  [4]. 
The  resonant  frequency  in  the  cavity  will  be  shifted  upward  in 
the  presence  of  plasma  and  the  measured  frequency  shifting  will 
determine  the  plasma  density  n.  Fig.  2  in  [4]  is  the  microwave 
interferometer  including  the  plasma  loaded  cavity.  Two  detected 
signals  //  =  |T|~  and  r  =  |T|sm  y  are  displayed  on  a  digital 
oscilloscope,  where  \T\  and  y  are,  respectively,  the  magnitude 
and  the  phase  of  the  transmission  coefficient  through  the  plasma 
cavity.  Assuming  uniform  plasma  in  the  cavity,  the  plasma  den¬ 
sity  n  is  expressed  by  [4] 
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Fig.  1.  Cross-sectional  view  of  X-band  TEn,t  mode  plasma  cavity.  A  thin 
stainless-steel  mesh  divides  the  cylindrical  plasma  container  into  the  plasma 
cavity  at  the  upper  portion  and  the  discharge  chamber  for  plasma  production  at 
the  lower  portion.  High-voltage  pulse  is  applied  at  the  terminal  at  the  bottom. 
Inset  shows  a  typical  current  waveform  for  pulsed  discharge  measured  by  a 
Rogowski  coil.  Distance  between  the  position  of  plasma  production  and  the 
antennas  is  2.8  cm. 


(a) 


where  F  is  the  sensing  microwave  angular  frequency  in  giga¬ 
hertz,  D(tf  2ir  is  the  half  half-width  of  the  resonant  curve  of  the 
cavity  without  plasma  in  megahertz,  and 


is  a  constant  with  the  mass  m  of  charged  carriers  which  respond 
to  the  microwave  fields.  In  the  case  of  gaseous  plasma,  m  is 
the  mass  of  an  electron  9.11  x  10~31  kfi,  while  the  mass  m  in 
(2)  may  be  given  by  the  effective  mass  of  a  p-ion  in  Llle,  or 
-.07  x  1()~2°  kg  [6]-[8],  40  limes  the  mass  of  a  4  IK*  atom  that 
is  smaller  than  the  mass  of  n-ions.  When  the  cavity  is  filled  with 
LHe,  the  charged  carriers  responding  to  the  microwave  fields 
may  be  p-ions  that  have  effective  mass  *2.9  x  103  times  that  of 
electrons,  resulting  in  the  constant .  1  =  7.29  x  101*  in  (2),  while 
.1  =  2.19  x  1013  for  the  single  electron  mass.  Using  (l)  and  (2) 
with  the  measured  signals  a  and  r  at  each  time  after  the  pulsed 
discharge,  one  can  obtain  effective  plasma  density  n  averaged 
over  the  cavity.  This  is  the  methodology  of  the  present  cavity 
measurement  of  plasma  in  Llle. 

Ill,  Experiment 

A  Experimental  Setup  and  Procedure 

The  schematic  cross  section  of  the  X-band  stainless-steel 
cylindrical  cavity  for  TEm  mode  with  resonant  frequency 


Fig.  2.  Measurements  in  cryogenic  helium  gas  with  temperature  1.69  K  and 
applied  voltage  less  than  20  kV.  (a)  Trace  of  microwave  signals  n  and  r  through 
the  cavity  when  it  is  tilled  with  cryogenic  helium  gas  with  pressure  0.3  torr. 
(b)  Logarithmic  plot  of  plasma  density  versus  time  when  the  cavity  is  filled 
witli  helium  gas  for  three  pressures  0.063,  0.070,  and  0.080  torr.  Inset  shows 
the  current  waveform  for  pulsed  discharge.  Decay  lime  constant  and  ambipolar 
diffusion  coefficient  are  obtained  from  the  fitted  straight  line. 

9.20  GHz  is  shown  in  Fig.  1.  The  sizes  of  the  cavity  are  3.6 
cm  inner  diameter  and  5.6  cm  in  length.  A  thin  and  flat  stain- 
less-steel  mesh  divides  the  cylindrical  vessel  into  two  portions: 
the  upper  for  the  plasma  cavity  of  microwave  measurement  and 
the  lower  for  the  discharge  chamber  for  plasma  production.  The 
mesh  separates  the  cavity  from  the  place  of  high-voltage  pulse 
application  to  avoid  the  damage  to  the  microwave  measurement 
circuit.  The  wall  and  two  antennas  in  the  cavity  are  grounded  to 
minimize  induction  noise  caused  by  the  high-voltage  pulse.  In 
the  discharge  chamber  with  an  inner  diameter  2.2  cm  and  the 
height  3.0  cm,  the  plasma  is  produced  by  single  or  repetitive 
shots  of  a  high-voltage  pulse  between  a  tungsten  needle  elec¬ 
trode  with  tip  curvature  50  /mi  and  ihe  stainless-steel  mesh.  The 
separaiion  between  the  needle  and  the  mesh  is  less  than  1  mm. 
The  pulse  voltage  can  be  varied  from  16  to  20  kV.  A  typical 
wave  form  of  a  discharge  current  measured  by  a  Rogowski  coil 
is  shown  in  the  inset  at  right  bottom  of  Fig.  1.  The  maximum 
discharge  current  is  less  than  360  A  and  the  duration  is  7  //s. 
While  (I)  is  applicable  only  to  a  uniform  plasma  as  in  gas 
helium,  we  expect  dial  the  plasma  density  in  LHe  averaged  over 
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the  cavity  volume  in  the  late  afterglow  period  may  be  evaluated 
reasonably,  if  the  plasma  frequency  is  less  than  the  microwave 
frequency  in  the  cavity. 

The  plasma  cavity  is  connected  to  a  microwave  circuit  by 
two  semi-rigid  cables  (UT-141A)  through  a  pair  of  antennas 
at  the  side  wall.  The  locations  and  shapes  of  the  antennas  are 
fixed  to  obtain  optimal  performance  for  the  measurements  of 
the  resonant  cavity  after  many  trials  and  error  tests.  Both  an¬ 
tennas  for  input  and  output  ports  are  roughly  identical  and  lo¬ 
cated  symmetrically.  A  CW  Gunn  oscillator  with  microwave 
power  less  than  1  mW  is  connected  to  the  input  antenna.  The 
output  antenna  is  connected  to  a  coaxial  detector  (HP432B)  and 
a  mixer  (Watkins-Johnson,  M80C),  as  shown  in  Fig.  2  in  [4]. 
The  cavity  is  installed  inside  a  LHe  dewar  of  volume  50  L  made 
of  stainless-steel  with  an  inner  diameter  45  cm  and  the  depth 
53  cm.  The  LHe  dewar  is  installed  inside  a  vacuum  container  in 
a  stainless-steel  cryostat  with  a  small-scale  refrigerator  to  pro¬ 
duce  LI  le  a  few  liters  per  hour.  The  whole  system  including  the 
TEnrj  cavity  is  evacuated  by  a  6-in  turbo-molecular  pump  to  a 
pressure  less  than  10"°  ton*  before  operation  of  the  refrigerator. 
Microwave  measurements  of  plasma  density  are  carried  out  for 
LHe  with  saturated  vapor  pressure  and  cryogenic  helium  gas  at 
a  temperature  range  from  T  =  4.2  to  1.69  K.  There  are  five 
small-radius  holes  at  the  top  wall  of  the  cavity  shown  in  Fig.  1. 
Through  the  holes,  one  can  watch  with  a  TV  camera  the  bright¬ 
ness  of  the  pulsed  discharge  through  the  viewing  port  at  the  top 
of  the  cryostat.  The  signal  from  the  Rogowski  coil  for  moni¬ 
toring  discharge  current  is  used  to  trigger  the  traces  on  the  dig¬ 
ital  oscilloscope  to  observe  the  waveforms  of  the  signals  u  and 
r  through  the  plasma  cavity.  In  our  experiment,  the  lower  limit 
of  temperature  of  1.69  K  is  attained  by  a  rotary  pump  (30  nr3/li) 
with  evacuation  cooling  of  the  50  L  LHe  dewar  surrounding  the 
plasma  cavity.  It  takes  8  h  to  cool  down  the  cavity  from  4.2  to 
1 .69  K,  and  around  50%  of  50  L  LHe  in  the  dewar  must  be  evap¬ 
orated  to  attain  2.0  K. 

The  resonant  frequency  is  measured  to  be  uv,/27r  = 

0.00-9.20  GHz  and  the  loaded  Q  value  of  the  TE11;i  mode 
cavity  without  plasma  is  determined  as  Q  =  u,y» /'2Bo  ~  1540, 
as  shown  in  Fig.  1  in  [4].  Because  of  the  slow  variation  in 
the  signals  u  and  v  on  the  oscilloscope,  the  amplitudes  of  the 
signals  and  setting  of  the  frequency  at  u.y,  -I-  B„  have  to  be 
calibrated  and  adjusted  before  each  shot, 

B,  Experimental  Results 

In  Section  ITI-B-L  the  measurements  of  decaying  plasma  in 
cryogenic  helium  gas  are  described  briefly,  and  the  results  on 
decaying  plasma  in  LHe  are  reported  in  Section  III-B-2. 

1)  Decaying  Plasma  in  Cryogenic  Helium  Gas  at 
T  —  4.2-7.69  K:  Fig.  2(a)  shows  waveforms  of  the  sig¬ 
nals  u  and  v  through  the  plasma  cavity  in  which  helium  gas  at 
a  temperature  1.69  K  and  pressure  0.3  torr  is  filled.  The  two 
signals,  0  <  u  <  1.0  and  —0.-5  <  v  <  0.5,  are  adjusted  to 
the  range  of  six  vertical  divisions,  respectively.  The  temporal 
evolution  of  the  logarithmic  plasma  density,  calculated  by  ( 1 ),  is 
shown  in  Fig.  2(b).  The  observed  plasma  density  is  in  the  range 
of  10*  -  10,J  cm’Jt  while  the  decaying  plasma  in  a  cryogenic 


helium  gas  in  a  TEnn  mode  S-band  cavity  [4]  was  character¬ 
ized  by  the  density  in  the  range  of  10°  -  10*  cm-3.  Ambipolar 
diffusion  was  identified  to  be  the  dominant  loss  mechanism  of 
the  decaying  plasma  in  the  S-band  TEnn  mode  cavity.  The 
present  TE113  X-band  cavity  has  a  sensing  frequency  about 
three  times  larger  and  Q  value  about  1/3  times  smaller  than  the 
S-band  TEon  mode  cavity.  Because  higher  plasma  densities 
are  observed  in  the  present  experiment,  we  expect  the  plasma 
to  decay  by  diffusion  as  well  as  by  recombination.  When  the 
pressure  is  less  than  0.085  torr,  the  loss  mechanism  is  dom¬ 
inated  by  diffusion  because  of  the  large  ambipolar  diffusion 
coefficient,  while  the  recombination  would  play  a  major  role 
in  the  loss  mechanism  for  the  pressures  higher  than  0.085  torr. 
The  temporal  plasma  decay  is  shown  in  Fig.  2(b)  for  pressures, 
0.063,  0.070,  and  0.080  torr  at  1.69  K,  where  the  current  wave 
form  is  depicted  in  the  inset.  We  note  that  the  peak  current, 
1.8  A,  is  much  less  than  that  shown  in  the  inset  of  Fig.  1  for 
the  discharge  in  LHe.  This  is  because  the  plasma  production 
in  the  cryogenic  gas  is  much  easier  than  in  LHe.  As  is  shown 
in  Fig.  2(b),  the  density  decreases  exponentially  with  Lime, 
suggesting  the  presence  of  the  ambipolar  diffusion.  The  decay 
time  increases  with  pressure,  as  is  expected  by  the  increase  in 
collision  frequency  between  helium  ions  and  neutral  helium 
atoms.  We  evaluate  the  ambipolar  diffusion  coefficient  Dn  from 
the  measured  decay  time  r  through  the  relation  r  =  A -/£>„, 
where  A  =  ((2.405//?)- +  (tt/I)2]-1/-  =  0.69nn  with  radius 
/?  =  1.8  cm  and  length  L  =  5.0  cm  for  the  cylindrical  cavity. 
The  diffusion  coefficients,  thus  obtained,  are  26.5,  11.3,  and 
8.66  onr/s  for  0.063,  0.070,  and  0.080  torr,  respectively.  The 
ambipolar  diffusion  coefficients  in  a  gas  with  pressure  P  may 
be  given  by  D A  =  540/P  (cm2/s)  for  the  helium  gas  at  300  K 
(see  Brown  [9]).  Assuming  D^P  =  540(7/300) i/:2,  with  a  gas 
temperature  T  in  equilibrium  T  =  7I>  =  T\  =  1.09  K  and  the 
momentum  transfer  collision  cross-section  5.8  x  10“ 1,7  enr  be¬ 
tween  helium  ions  and  neutral  helium  atoms,  we  find  expected 
values  of  Dn  to  be  3.6,  3.3,  and  2.9  enr/s  at  pressures  0.063, 
0.070,  and  0.080  torr,  respectively.  Although  the  observed 
values  of  are  somewhat  larger  than  the  expected  values,  a 
deviation  might  be  caused  by  an  additional  loss  of  the  plasma 
such  as  the  expansion  of  the  plasma  through  the  stainless  mesh 
from  the  cavity  to  the  discharge  chamber  where  the  plasma  was 
originally  produced.  The  late  afterglow  plasma  in  the  present 
cavity  disappears  faster  than  a  case  for  plasma  cavity  enclosed 
completely  with  metal  wall,  as  was  discussed  in  [4]. 

Fig.  3(a)  shows  the  temporal  decay  of  plasma  density  at  T  - 
LG9  K  in  helium  gas  for  pressures  larger  than  those  shown  in 
Fig.  2(b).  The  decay  is  slower  for  the  case  of  0.090  torr  than  the 
case  of  0.220  torr.  Such  a  pressure  dependence  does  not  agree 
with  the  diffusion  process  described  by  D^P  —  const.,  and 
we  see  that  the  decay  lines  are  not  straight  lines  in  the  semi- log 
plot,  as  shown  in  Fig.  3(a).  We  may  conclude  that  the  loss  of  the 
plasma  is  caused  not  by  ambipolar  diffusion,  but  by  recombina¬ 
tion  [9],  [10].  Fig.  3(b)  shows  the  inverse  of  the  plasma  density 
in  Fig.  3(a)  versus  time.  The  slope  of  the  fitted  straight  lines 
gives  the  values  of  recombination  coefficient  n.  Estimated  n 
thus  obtained  are  1.9  x  10“7, 2.8  x  10“7  and  7.0  x  10“7nn3/s 
for  pressure  0.090,  0.145  and  0.220  torr,  respectively.  The  ob¬ 
served  recombination  coefficient  #v  decreases  with  temperature 
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Fig.  3.  Temporal  evolution  of  the  plasma  density  when  the  cavity  is  filled  with 
helium  gas  at  temperature  1 .69  K  for  three  gas  pressures  0.090, 0. 145,  and  0.220 
torr.  Applied  voltage  is  less  than  20  kV.  (a)  Logarithmic  plasma  density  versus 
time,  (b)  Inverse  of  plasma  deasity  versus  time.  Recombination  coefficient  is 
obtained  from  the  fitted  straight  line. 


T,  agreeing  with  the  previous  work  [  10],  and  increases  with  gas 
pressure  as  expected  for  the  ion-ion  recombination  coefficient 
given  by  Thomson’s  theory  for  tenuous  gases  [11]. 

2)  Decaying  Plasma  in  the  Cavity  Filled  With  Liquid  Helium 
at  T  =  4.2 -1.69  K:  Next,  the  cavity  is  filled  with  LHe,  and  the 
signals  u  and  r  are  measured.  The  experiments  are  carried  out 
in  LHe  mostly  at  saturated  vapor  pressure  under  the  condition 
that  the  gas  volume  above  the  free  surface  of  LHe  filling  conr 
pletely  the  cavity  is  isolated  from  the  outer  helium-gas-feeding 
tubing-connection  by  closing  the  valve.  Since  LHe  in  the  cavity 
is  in  thermal  equilibrium  with  the  surrounding  LHe  in  the  dewar 
that  is  cooled  down  by  evacuation,  LHe  in  the  cavity  is  main¬ 
tained  to  have  the  saturated  vapor  pressure.  The  temperature  of 
LHe  in  the  cavity  is  measured  precisely  from  the  gas  pressure 
above  the  free  surface.  Plasma  density  is  calculated  by  (1)  in 
which  bare  electrons  are  assumed  to  respond  to  the  microwave 
fields.  The  signals  //  and  r  on  the  order  of  200  ms  after  the 
pulsed  discharge  are  observed  with  good  reproducibility,  while 
the  plasma  decays  in  much  shorter  period  in  the  gas  case  in 
Figs.  2  and  3.  The  long  time  sustainment  of  the  plasma  may 
result  from  the  slow  diffusion  of  the  localized  plasma  through 
LHe  compared  to  the  gas  cases  given  in  the  previous  sub-sec¬ 
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Fig.  4.  Measured  responses  through  the  plasma  cavity  filled  with  liquid  helium 
at  saturated  vapor  pressure  for  various  temperatures  with  applied  voltage  20  kV. 
(a)  Waveforms  of  signals  «  and  r.  (b)  Calculated  plasma  density  versus  time  for 
two  temperatures  2.87  and  3.49  K. 

tion.  Fig.  4(a)  shows  the  responses  of  the  signals  u  and  e  through 
the  plasma  cavity  for  various  temperatures  with  the  applied  dis¬ 
charge  voltage  of  20  kV.  The  plasma  density  is  estimated  to  be 
on  the  order  of  10*  mi”\  an  approximate  value  in  the  late  af¬ 
terglow  period  of  the  pulsed  discharge  in  LHe.  Temporal  evo¬ 
lutions  of  the  plasma  density  calculated  by  (1)  are  shown  in 
Fig.  4(b)  for  temperatures  2.87  and  3.49  K.  The  maximum  den¬ 
sity  is  attained  at  T  -  2.87  K,  while  the  plasma  signals  disap¬ 
pear  below  T  —  2.3  K.  The  dependence  of  the  u  and  r  signals 
for  various  applied  voltages  in  the  case  T  =  2.87  K  is  shown  in 
Fig.  5(a),  and  the  corresponding  plasma  densities  are  shown  in 
Fig.  5(b)  for  applied  voltages  17,  19  and  20  kV.  It  is  shown  that 
the  plasma  density  increases  with  the  applied  voltage.  The  range 
of  observed  plasma  density  is  !<>'  -  10iJ  mi”3  in  the  cavity, 
which  is  estimated  with  the  assumption  that  the  charged  car¬ 
riers  responding  to  the  microwave  fields  in  the  cavity  are  single 
electrons  rather  than  ion  clusters  [6]-[8]. 
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Fig.  5.  Measured  responses  through  the  plasma  cavity  filled  with  liquid  helium 
at  saturated  vapor  pressure  for  various  applied  voltages.  Liquid  temperature 
is  2.87  K  at  saturated  vapor  pressure,  (a)  Waveforms  of  signals  it  and  r.  (b) 
Calculated  plasma  density  versus  time  for  three  applied  voltages  17,  19,  and  20 
kV. 

IV.  Analysis  of  Temporal  Evolution  of  Localized 
Plasma  in  Liquid  Helium 

Simple  three-dimensional  analysis  of  localized  decaying 
plasma  in  LHe  is  presented  in  this  section.  Recombination  co¬ 
efficients  and  mobility  of  ions  in  LHe  available  in  the  Literature 
are  reviewed  and  expected  temporal  evolution  of  the  decaying 
plasma  in  LHe  is  analyzed  numerically. 

A.  Properties  of  Charged  Carriers  in  Liquid  Helium 

Properties  of  p-  and  n-ions  in  LHe  have  been  studied  by  many 
researchers  [12]t  revealing:  )  the  p-ion  in  LHe  to  be  a  semi- 
macroscopic  impurity  object  with  radius  of  0.6  nm  and  with  ef¬ 
fective  mass  40  times  that  of  a  4IIe  atom  for  the  temperature 
range  1. 2-4.2  K  and  2)  the  n-ion  in  LHe  to  be  a  semi-macro¬ 
scopic  impurity  object  with  radius  1-1.6  nm  and  with  effective 
mass  240  times  that  of  a  4Ho  atom  [6]-[8L  [13].  Ambient  neu¬ 
tral  4Ho  atoms  are  attracted  by  a  week  dipole  field  of  a  core  4H<- 
p-ion,  forming  a  p-ion  cluster  or  an  ice-ball.  On  the  other  hand, 
a  free  electron  in  LHe  is  trapped  inside  a  vacuum  bubble  which 
is  explained  by  Pauli’s  exclusion  principle.  Because  of  the  hy¬ 
drodynamic  nature,  the  electron  in  the  bubble  is  called  the  n-ion. 
At  0.3  K  and  24  atm,  p-  and  n-ions  can  move  without  drag  force 
when  the  velocities  are  below  the  Landau  limit  46  m/s  known  to 
create  a  pair  of  rotons  [14].  Large  ion  mobility  was  observed  in  a 
variety  of  experiments  in  LHe  for  higher  temperatures  and  lower 
pressures  with  the  ion  velocity  much  smaller  than  the  Landau 


velocity.  So  far,  all  of  the  experimental  studies  on  ions  in  LHe 
have  been  restricted  to  direct  current  (dc)  current  densities  much 
less  than  10  “r‘  A/cm"  produced  by  a  radio-isotope  ion  source 
16]-[8],  [15],  [16]  or  dc  field  emission  from  a  metal  projection 
1 14]  immersed  in  LHe.  The  masses  Mv  and  A/„  of  p-  and  n-ions 
are  roughly  M},  =  40 M(4H<  )  =  2.08  x  I0“2"  kg  and  M,t  = 
240 Mi 4 Ilf  )  =  1.01  x  10--4  kg, respectively  [6]-[8].  Since  the 
mass  of  a  p-ion  is  smaller  than  that  of  the  n-ion,  relatively  mo¬ 
bile  charged  carriers  that  respond  to  the  microwave  field  in  our 
experiment  may  be  mainly  p-ions.  In  LHe,  the  responding  p-ion 
has  the  mass  Mv  that  is  MP/M  =  10  x  4  x  1837  =  2.0  x  10*’ 
times  the  mass  Mc  of  a  single  electron.  The  range  of  observed 
plasma  density  may  be  2.9  x  10°  times  larger  than  those  in  the 
plasma  experiment  in  cryogenic  helium  gas  shown  in  the  pre¬ 
vious  subsection,  because  the  observed  plasma  density  n  is  pro¬ 
portional  to  mass  A/p  for  a  given  plasma  frequency. 

The  ion-ion  recombination  coefficient  n  in  LHe  was  mea¬ 
sured  by  Careri  and  Gaeta  [15]  for  a  temperature  range  0.87-2.0 
K  at  saturated  vapor  pressure,  as  shown  by  white  squares  in 
Fig.  6(a)  where  theoretical  values  are  also  shown.  The  measured 
a  was  found  to  agree  with  the  conventional  Langevin’s  theory 
shown  by  black  circles  for  high  pressure  gases  at  300  K,  where 
the  mean  free  paths  of  ions  are  much  smaller  than  the  average 
ion  inter-particle  distance.  Thus,  the  recombination  coefficient 
may  be  given  by  [9]  n  =  t  (pr  +  //„)/, r«,,  where  pp  and  //„ 
are  the  mobilities  of  the  p-  and  the  n-ions,  and  o(em3/s)  = 
1.8ixi0“c(//p+y/„  )  with//p  and y/„  intheunitof  cm-s_1V"1. 
Using  the  data  of  the  mobility  by  Schwarz,  the  black  circles 
of  a  in  Fig.  6(a)  denoted  by  Langevin-Schwarz  are  obtained. 
White  triangles  in  Fig.  6(a)  are  Harper's  theory  [17]  of  a  in 
high  pressure  gases  obtained  in  somewhat  improved  way  rather 
chan  Langevin's  theory  and  the  results  are  around  1/4  times 
smaller  than  those  of  Langevin-Schwarz.  The  white  triangles 
are  denoted  as  Harper-Schwarz,  since  the  data  by  Schwarz  were 
again  used.  The  values  of  o  measured  by  Careri  and  Gaeta  [15] 
are  located  just  between  two  curves  of  Langevin-Schwarz  and 
Harper-Schwarz.  The  mobility  of  p-ions  in  LHe  is  determined 
by  the  interactions  of  charged  carriers  with  rotons,  phonons  and 
Hlo  impurities  that  are  the  targets  of  scattering  for  the  collisions 
[18].  These  exitons  play  the  role  of  neutral  atoms  to  provide  the 
mobility  of  ions  in  conventional  gases.  For  T  above  0.6  K,  the 
dominant  exitons  in  LHe  are  rotons,  and  Langevin’s  recombi¬ 
nation  coefficient  is  given  by 

a  (cm'3  /s)  =  Ux  10  “  J  oxp(  8.7 /T )  (3) 

that  is  denoted  as  “theory  of  roton"  by  a  broken  line  in  Fig.  6(a). 
The  measured  recombination  coefficient  by  Careri  and  Gaeta 
[15]  agrees  well  with  the  broken  line  of  roton  for  T  <  2.0  K. 
We  extend  the  recombination  coefficient  n  given  by  (3)  to  the 
temperature  range  T  =  0.-3-2.17  K,  although  no  experimental 
data  measurements  were  available  in  the  full  range.  The  white 
circles  at  the  left  top  comer  of  Fig.  6(a)  are  ion-ion  recombi¬ 
nation  coefficients  n  given  by  Thomson’s  theory  [11]  for  low 
gas  pressure  less  than  1  atm  such  that  mean  free  paths  of  ions 
are  much  larger  than  the  average  ion  inter-particle  distance.  The 
ion-ion  recombination  coefficient  is  n  ^  10“-  cm4/s  around 
temperature  T  =  0.5  K  as  shown  in  Fig.  6(a). 


6 


IEEE  TRANSACTIONS  ON  PLASMA  SCIENCE 


(a) 


Fig.  6.  Recombination  and  diffusion  coefficients  of  ions  in  liquid  helium 
calculated  from  known  data  in  the  literature,  (a)  Recombination  coefficients 
versus  temperature  measured  by  Careri  and  Gaeta  [15].  Other  curves  arc 
theoretical  values  using  the  data  of  mobility  of  positive  and  negative  ions 
measured  by  Schwarz  [16].  (b)  Free  and  ambipolar  diffusion  coefficients  versus 
temperature  calculated  from  the  data  of  mobility  measured  by  Schwarz  [16]. 


Schwarz  measured  zero-field  mobility  of  tenuous  p-  and 
n-ions  in  LHe  at  saturated  vapor  pressure  [16].  The  mobility 
of  the  p-ion  at  T  =  0.37  K  increases  as  much  as  10°  times 
the  lowest  values  at  4.2  K,  while  the  mobility  of  the  n-ion  at 
T  =  0.37  K  increases  as  much  as  3  x  104  times  the  lowest 
values  at  4.2  K.  The  ambipolar  diffusion  coefficient  D a  of 
ions  in  LHe  may  be  evaluated  by  Dp.»  =  or 

D;i  ff(cnr/s)  =  8.63  x  lO”*>T(A')/iP.ll(cm“/Vs),  Here,  we 
assume  Tv  =  Tn  =  T.  The  free  diffusion  coefficients  as  a 
function  of  temperature  are  shown  by  broken  lines  in  Fig.  6(b). 
The  ambipolar  diffusion  coefficient  An  may  be  given  by  [9] 


A.  (onr/s) 


_  Pn&p  “1"  Pp&n 

/'»  +  fip 

=  1.73  x  l()-4T(K)ftil(cm2/Vs) 


(4) 


where  i  jftu  -  i/p„  +  1  j\iv  as  given  by  Blanc's  Law.  The 
ambipolar  diffusion  coefficient  An  is  shown  by  a  solid  curve  in 
Fig.  6(b).  It  is  evident  that  Dn  <  An  <  Dv,  although  these 
three  diffusion  coefficients  are  on  the  same  order  of  magnitude 
for  T  >  1.0  K.  As  shown  in  Figs.  6(a)  and  6(b),  both  An  and  n 
increase  many  orders  of  magnitude  wilh  the  decrease  in  T  from 
2  to  0.5  K  because  of  the  nature  of  superfluidity. 


B.  Evolution  of  Localized  Plasma  in  Liquid  Helium 

The  localized  plasma  produced  by  a  pulsed  discharge  in  LHe 
is  characterized  by  the  high  initial  density  and  the  tempera¬ 
ture  on  the  order  of  101*  cm_;j  and  3  x  104  K,  respectively, 
surrounded  by  a  gas  layer  [2],  [5].  When  the  pulsed  discharge 
is  over,  the  plasma  vanishes  quickly,  and  some  ions  are  im¬ 
planted  into  LHe  before  they  recombine  [5].  Ions  with  tempera¬ 
ture  higher  than  the  temperature  of  LHe  will  diffuse  quickly  and 
may  distribute  through  LHe  in  the  cavity. 

Decay  in  density  n{r.  /  )  of  the  localized  plasma  produced  in 
LHe  may  be  caused  by  recombination  and  diffusion  losses  as 
expressed  by 
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where  n  and  Du  are  the  recombination  coefficient  and  the  am¬ 
bipolar  diffusion  coefficient,  respectively.  The  differential  equa¬ 
tion  (5)  is  solved  numerically  assuming  the  initial  density  distri¬ 
bution  observed  earlier  in  Fig.  3(a)  in  [5],  namely,  the  distribu¬ 
tion  wilh  the  half  half- width  of  0.3  mm.  First,  we  choose  the  co¬ 
efficients  n  and  A,  from  Fig.  6  and  set  o  =  1.2  x  10“ f  cm3/s 
and  An  =  L0  x  lll”°  ex lr/s  in  a  normal  fluid  case  ( Tp  = 
T„  -  T  —  1.2  K),  while  n  =  1.0  x  10“*'  cm3/s,  Da  =  1.26  x 
10“4  cu r/s  in  a  superfluid  case  (Tv  =  Tu  =  T  =  1.0  K).  The 
results  for  both  cases  are  shown  in  Fig.  7(a)  and  (b),  respectively. 
As  shown  in  Fig.  7(a).  the  effect  of  ambipolar  diffusion  becomes 
dominant  for  the  density  less  than  10"  cm-"  at  the  time  100  s  in 
the  normal  fluid,  and  1M0  s  in  the  superfluid  case  in  Fig.  7(b). 
In  our  experiment,  as  shown  in  Figs.  4(b)  and  5(b),  however,  the 
plasma  density  decays  in  less  than  1  s.  Such  a  disagreement  be¬ 
tween  observed  and  expected  decay  times  may  be  attributed  to  a 
difference  in  the  diffusion  coefficients  between  the  experiment 
in  Figs.  4(b)  and  5(b),  and  the  calculation  in  Fig.  7. 

Next,  we  choose  values  of  the  diffusion  coefficients  from  the 
following  expression  as 


A, 


2  371, 


(6) 


where  L  —  3.6  cm  is  the  length  of  the  cylindrical  plasma  con¬ 
tainer  and  r.i  =  0.2  s  is  the  delay  time  after  the  pulsed  dis¬ 
charge.  The  decay  time  T<\  in  (6)  is  the  time  of  signal  arrival  at 
the  location  of  the  antenna  at  r  —  Lj'l  from  a  point  source  of 
plasma  at  r  =  0.  The  diffusion  coefficient  A*  —  6.3  cur/s 
obtained  from  (6)  is  found  to  agree  with  the  observation  in 
the  decaying  plasma  in  cryogenic  gas  shown  in  Fig.  8(a)  of 
|4],  Fig.  8(a)  in  this  paper  shows  the  temporal  evolution  of  the 
plasma  density  versus  the  radial  distance  r  assuming  D a  = 
G.5  cni'/s  and  n  =  1.0  x  10“ 1  cm^/s.  We  note  that  plasma  ex¬ 
pands  much  faster  than  the  expected  results  shown  in  Figs.  7(a) 
and  (b).  To  see  the  effect  of  the  diffusion  coefficients  on  the 
plasma  density  at  the  antenna,  we  plot  the  plasma  density  at  the 
antenna  location  r  =  1/2  =  2.8  cm,  as  a  function  of  lime  by 
the  solution  of  (5)  with  n  =  ft 
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Fig.  7.  Temporal  evolution  of  localized  cryogenic  plasma  in  liquid  helium 
with  a  diffusion  coefficient  calculated  by  (5).  Here,  an  initial  density  profile 
measured  by  Kojima  el  ai  [5]  and  recombination  and  diffusion  coefficients 
shown,  respectively,  in  Fig.  6(a)  and  (b)  have  been  used,  (a)  Liquid  helium  at 
constant  temperature  T  =  Tr  =  Tts  =  4.2  Iv  (normal  fluid)  with  saturated 
vapor  pressure,  (b)  Liquid  helium  at  constant  temperature  T  —  Tt.  =  Tn  = 
1.0  K  (superfluid)  with  saturated  vapor  pressure.  Here  T,  Tr,  and  T„  are  the 
temperature  of  liquid  helium,  p-ions,  and  n-ions,  respectively. 

where  C  is  a  constant.  The  origin  r  =  0  is  taken  at  the  elec¬ 
trode  where  the  plasma  is  produced.  Fig.  8(b)  shows  the  plasma 
density  as  a  function  of  time  at  the  location  r  =  2.8  <  m  for 
Ai  =  2.9  ~  G.o  enr/s.  The  density  becomes  maximum  at 
the  time  when  On/Oi  =  0  in  (7)  that  results  in  (6),  which  cor¬ 
responds  to  the  delay  time  T<\.  Pig.  8(b)  would  compare  favor¬ 
ably  with  the  experimental  data,  Pig.  4(b)  and  Fig.  5(b),  in  the 
time  scale  except  the  observed  oscillatory  nature  with  three  or 
four  peaks.  The  observed  diffusion  process  takes  place  much 
faster  than  those  expected  from  Fig.  7,  where  heavy  p-ion  mass 
was  assumed.  The  temporal  observation  of  the  plasma  density 
as  shown  in  Figs.  4(b)  and  5(b)  might  be  understood  by  the  loss 
process  with  the  diffusion  coefficient  given  by  (6),  rather  than 
by  (4).  We  might  conclude  that  the  ambipolar  diffusion  involves 
ions  with  mass  of  light  single  helium  atom  and  bare  electrons  in 
LHe.  The  oscillatory  nature  of  the  plasma  density  observed  in 
Pig.  4(b),  as  well  as  in  Fig.  5(b)  with  a  number  of  peaks  remains 
to  be  solved. 

V.  Discussion  and  Conclusion 

We  have  measured  the  density  of  the  decaying  plasma  in 
LHe.  The  observed  masses  of  charged  carriers  are  estimated 
to  be  much  smaller  than  those  expected  in  [6J— [8].  We  con¬ 
clude  that  the  recombination  coefficient  [15]  and  the  mobility 


fa) 


Fig.  8.  Numerical  solutions  for  the  evolution  of  a  localized  plasma  in  liquid 
helium  with  assumed  diffusion  coefficients  to  lit  the  observed  results  in 
Figs.  4(b)  and  5(b).  (a)  Temporal  evolution  of  localized  cryogenic  plasma  in 
liquid  helium  with  initial  density  distribution  measured  by  Kojima  el  ai  [5). 
Here,  the  recombination  coefficient  =  1.0  x  lfl-7  cm '/s  and  ambiploar 
diffusion  coefficient  Dn  =  G.o  cm:/s  are  assumed,  (b)  Temporal  evolution 
of  plasma  density,  (7),  at  the  position  of  the  antennas  r  =  1/2  =  2.8  rm 
from  the  point  source  of  plasma  discharge  at  r  =  0  and  t  =0.  Recombination 
coefficient  =  0  and  ambiploar  diffusion  coefficients,  D;.  =  G.o,  4.4  and 
2.9  cm2/*  arc  assumed.  This  figure  is  to  be  compared  with  those  observed  in 
Figs.  4(b)  and  5(b). 

[16]  measured  previously  at  steady-state  in  drift  tubes  filled 
with  LHe  using  dc  ion  source  may  not  be  applied  to  the  tran¬ 
sient  state  just  after  the  pulsed  discharge  in  the  present  exper¬ 
iment.  Our  microwave  measurement  is  made  using  the  prin¬ 
ciple  of  comparing  the  sensing  frequency  u.*  with  plasma  fre¬ 
quency  u/r  =  {nr2 /com)L/2.  For  a  given  «; P ,  the  measured 
density  n  increases  with  mass  tn  in  (2).  Assuming  p-ions  to 
be  the  responding  charged  carriers  to  the  microwave  field,  the 
scales  of  density  n  in  vertical  axes  in  Figs.  4(b)  and  5(b)  would 
be  2.9  x  IIP  times  larger  than  those  shown  in  these  figures. 
Although  the  exact  masses  of  ions  in  LHe  just  after  the  high- 
voltage  pulsed  discharges  have  not  been  reported  yet  in  any  lit¬ 
erature,  it  may  be  reasonable  to  estimate  the  density  of  charged 
carriers  to  be  on  the  order  of  10*  fin-*5  with  the  assumption 
of  a  single  electron  mass  in  the  late  afterglow  period  rather 
than  1013  cm-3  assuming  heavy  p-ions.  The  present  experi¬ 
ment  suggests  that  masses  of  charged  carriers  observed  after 
the  high-voltage  pulsed  discharge  in  LI  le  are  much  smaller  than 
those  of  heavy  ions  in  steady  state  measured  in  LHe  [6]— 18].  It 
is  conjectured  that  some  finite  time  on  the  order  of  one  second 
will  need  to  form  positive  ion  clusters  or  electron  bubbles  with 
full  effective  masses  [6] — [8].  All  the  earlier  experiments  to  mea¬ 
sure  the  ion  masses  were  made  for  charged  carriers  produced 
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by  the  steady-state  ion  sources.  In  other  words,  ions  were  cre¬ 
ated  a  long  time  before  the  measurements.  Our  measurement, 
however,  is  made  just  after  the  production  of  charged  carriers  in 
LHe.  Our  result  does  not  exclude  the  possibility  of  the  presence 
of  heavy  ions  that  could  not  be  detected.  The  lower  limit  of  tem¬ 
perature  of  LHe  in  our  experiment  is  1.69  K.  Fig.  6(a)  and  (b) 
suggests  that  distinct  features  of  cryogenic  ions  in  LHe  would 
appear  around  T  =  0.5  K. 

The  present  measurements  of  cryogenic  afterglow  plasma  in 
LHe  are  carried  out  by  a  microwave  cavity  method.  Additional 
measurements  for  localized  low  density  plasma  either  by  elec¬ 
trostatic  probe  or  by  wave  propagation  would  be  desirable  to 
confirm  the  results  in  the  present  microwave  measurements  that 
lack  spatial  resolution.  The  Langmuir  probe  was  used  to  deter¬ 
mine  the  plasma  characteristics  in  cryogenic  afterglow  plasma 
in  LHe,  giving  a  preliminary  result  of  unexpected  high  temper¬ 
ature  of  charged  carriers.  Further  efforts  to  improve  accuracy  of 
the  measurements  are  under  way. 

In  conclusion,  charged  carriers  with  masses  smaller  than 
the  expected  conventional  values  are  observed  in  the  transient 
plasma  in  LHe.  Although  the  microwave  cavity  method  may 
not  detect  the  localized  plasma  precisely,  our  measurements  to 
estimate  the  masses  of  ions  in  the  cryogenic  decaying  plasma 
in  LHe  are  reported  for  the  first  time  in  the  literature.  Study 
of  cryogenic  plasma  in  LHe  including  the  time  of  formation 
of  an  ion  cluster  and  an  electron  bubble  with  various  effective 
masses  is  a  new  subject  of  research  in  plasma  science.  Distinct 
effects  of  superfluidity  in  cryogenic  plasma  in  LHe  near  0.5  K 
are  to  be  clarified  experimentally  in  the  future.  Theoretical  and 
computational  analyzes  to  elucidate  the  novel  features  of  the 
cryogenic  plasma  in  LHe  should  be  promoted. 
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1.  Introduction 

We  produced  cryogenic  decaying  plasma  in  liquid  helium  (LHe)  as  well  as  in  cryogenic  gas 
helium  with  temperature  below  4.2  Kelvin  by  means  of  a  pulsed  discharge  [1],  A  plasma  was 
produced  in  cryogenic  helium  gas  in  a  TEon  mode  S-band  cavity  with  very  long  decay  time  on 
the  order  of  1  sec  [1],  while  high-voltage  pulsed  discharge  with  needle  electrodes  in  LHe 
produced  the  localized  plasma  with  density  more  than  1018cm'3  and  electron  and  ion  temperatures 
of  a  few  eV  [2].  Upon  termination  of  the  discharge,  the  plasma  remained  as  long  as  on  the  order 
of  10 pstc  in  LHe.  Some  of  ions  and  electrons  remained  longer  in  LHe  before  they  recombine. 
Such  a  cryogenic  plasma,  or  positive  and  negative  ion  mixture  dissolved  in  LHe,  is  studied  in 
detail. 


2.  Experiment 

We  apply  the  microwave  cavity  method  to 
measure  decaying  plasma  in  LHe  [3].  The 
schematic  cross  section  of  plasma  cavity  is 
shown  in  Figure  1.  Plasma  is  produced  in 
plasma  cavity  by  means  of  a  pulse  discharge. 
Late  afterglow  plasma  in  an  X-band  cavity 
filled  with  LHe  in  the  temperature  range  from 
4.2  to  1.69  K  is  studied.  Our  microwave 
measurement  is  based  on  the  principle  of 
comparing  the  applied  frequency  ruwith  the 
cavity  resonance  frequency.  The  resonant 
frequency  in  the  cavity  will  be  shifted  upward 
in  the  presence  of  plasma  and  the  detected 
frequency  shifting  will  determine  the  plasma 
density  n  through  plasma  frequency 
cop  -  (ne2 1  £0m)112 . 
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Figure  1.  Cross  section  of  plasma  cavity 


3.  Results  and  Conclusion 

A  slow  decaying  plasma  in  LHe  at  a  saturating  vapor  pressure  with  temperature  above  2.3  K  is 
detected  as  long  as  200  ms  after  the  discharge.  Plasma  density  vs.  time  for  two  temperatures  2.87 
and  3.49K  is  shown  in  Figure  2.  The  exact  masses  of  ions  in  LHe  just  after  the  high-voltage 
pulsed  discharges  have  not  been  reported  yet  in  any  literature.  It  is  estimated  that  the  reasonable 
density  of  charged  carriers  to  be  on  the  order  of  108  cm'3  with  the  assumption  of  a  single  electron 
mass  in  the  late  afterglow  period  rather  than  1013  cm'3  assuming  a  positive  ion.  The  present 
experiment  suggests  that  masses  of  charged  carriers  observed  after  the  high-voltage  pulsed 
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discharge  in  LHe  are  much  smaller  than  those  of  ions  in  a  steady  state  measured  in  LHe  [4,5], 

We  conclude  that  masses  of  charged  carriers  measured  previously  at  steady-state  in  drift  tubes 
filled  with  LHe  using  DC  ion  source  may  not  be  applicable  to  the  transient  state  just  after  the 
pulsed  discharge  in  the  present  experiment  It  is  conjectured  that  some  finite  time  on  the  order  of 
one  second  will  need  to  form  positive  ion  clusters  or  electron  bubbles  with  full  effective  masses 
[4,5].  We  note  that  all  the  earlier  experiments  to  measure  the  ion  masses  were  made  for  charged 
carriers  produced  by  the  steady-state  ion  sources.  Our  measurement  however,  is  made  just  after 
the  production  of  charged  carriers  in  LHe.  Our  result  does  not  exclude  the  possibility  of  the 
presence  of  heavy  ions  that  could  not  be  detected.  Our  estimated  density  based  on  the  cavity 
method  suggests  an  effective  mass  of  ions  much  different  from  that  expected  by  the  conventional 
theory.  A  simple  diffusion  loss  model  for  a  temporal  evolution  of  the  localized  plasma  produced 
by  the  pulsed  discharge  in  LHe  reveals  the  presence  of  charged  carriers  with  masses  much  lighter 
than  those  expected  in  LHe[3]. 


Figure  2.  Piasma  density  vs.  time  for  two  temperatures  2.87  and  3.49  K 
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Observation  of  dust  particles  in  a  plasma  under  cryogenic  environment 

M.  Shindo,  C.  Kojima,  0.  Ishihara 
Yokohama  National  University 

A  complex  plasma  in  a  cryogenic  environment  is  studied  experimentally.  A  plasma  is 
produced  in  a  vapor  of  liquid  helium  by  applying  the  low  frequency  AC  voltage  between  plane 
electrodes.  Fine  acrylic  particles  are  injected  into  a  plasma.  Fine  particles  are  negatively  charged 
because  of  the  higher  mobility  of  electrons  than  ions  hitting  on  the  surface  of  the  particles. 
Electrons,  ions  and  fine  particles  form  a  complex  plasma  where  the  overall  charge  neutrality  is 
satisfied.  The  complex  plasma  is  quite  different  from  the  ordinary  plasma  in  a  sense  that  the 
coupling  constant  of  fine  particles,  defined  by 
the  ratio  between  the  Coulomb  interaction 
energy  and  the  thermal  energy  of  fine  particles, 
is  much  larger  than  unity  and  fine  particles  are 
strongly  coupled  through  the  Coulomb 
interaction. 

Figure  1  shows  the  experimental  apparatus  in 
which  Liquid  helium  is  in  a  double  glass 
Dewar  of  10  cm  in  diameter  and  lm  in  length, 
and  the  liquid  temperature  is  cooled  down  to 
IK  by  evaporative  cooling.  The  pressure  at  the 
upper  part  of  the  Dewar  was  about  300  Pa. 

The  discharge  electrodes  of  5mm  apart  are 
placed  about  20cm  above  the  surface  of  liquid  01  . 

tig.l  Schematic  diagram  oi  dust  plasma  apparatus 

helium,  while  a  mesh  electrode  and  a  plate 

in  cryogenic  condition. 

electrode  of  5  cm  in  diameter  are  placed 

below  the  discharge  electrodes.  The  10kHz  AC  voltage  of  lkV  was  applied  to  the  electrodes. 
The  discharge  electrodes  are  protected  by  a  cylindrical  insulator  to  suppress  the  rapid  rise  of  the 
temperature  of  both  Liquid  helium  and  the  vapor  gas,  while  all  the  electrodes  are  in  an  acrylic  tube 
of  5cm  in  diameter  to  confine  the  particles.  The  plasma  diffuses  toward  the  surface  of  Liquid 
helium  through  a  5mm  hole  in  the  lower  electrode.  Electron  temperature  decreases  gradually 
toward  the  surface  of  Liquid  helium.  The  acrylic  particles  of  5  pm  in  diameter  were  dropped 
through  the  5mm  holes  at  the  center  of  the  upper  electrode  into  the  plasma.  The  He-Ne  laser  light 
irradiated  though  a  1cm- width  slit  on  the  silver  plated  glass  Dewar  confirmed  the  particles 
levitating  1cm  above  the  plate  electrode. 

This  work  is  supported  by  Asian  Office  of  Aerospace  Research  and  Development. 
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Plasma  production  in  a  cryogenic  condition 

C.  Kojima,  M.  Kugue,  T.  Maezawa,  M.  Shindo  and  O.  Ishihara 
Yokohama  National  University 

A  steady  production  of  a  plasma  in  a  cryogenic  environment  is  studied  experimentally.  After 
a  successful  production  of  a  plasma  in  a  gas  at  the  liquid  helium  temperature  (below  4.2K)  by  a 
pulse  discharge  (20kV  with  a  pulse  width  7//s)[l],  our  efforts  are  concentrated  on  producing  a 
cold  plasma  for  an  extended  period  of  time.  Our  attempt  to  produce  a  cryogenic  plasma  in  a 
steady  state  is  achieved  by  an  electrical  discharge  (AClOkHz,  ~lkV)  in  a  vapor  (300  Pa)  above 
the  liquid  helium  surface. 

The  electrical  discharge  in  the  vapor 
necessarily  brings  a  heat  in  the  environment 
resulting  in  the  boiling  of  the  liquid  helium. 

To  overcome  the  rise  of  the  surrounding 
plasma  temperature  as  well  as  the  instability 
of  the  surface  of  the  liquid  helium,  a  special 
purpose  discharge  unit  was  developed  as 
shown  in  Fig.  1  in  which  a  region  of  plasma 
production  is  isolated  from  the 
surroundings  by  a  small  glass  Dewar[2], 


Two  identical  plane  discharge  electrodes  Ay  A; 

of  22  mm  in  diameter  and  5mm  in  thickness 
have  coaxial  holes  at  the  center  through 

which  the  helium  gas  was  fed  from  the  n  j  Schematic  diagram  of  the  discharge 
upper  hole.  A  plasma  is  produced  by  a 

discharge  between  the  electrodes  and  the  unit  m  the  vaPor  of  llciuld  hellum- 
gas  flow  takes  the  plasma  outside  of  the 
discharge  region  through  the  bottom  hole. 

The  rise  of  the  plasma  temperature  once  outside  of  the  glass  Dewar  is  successfully  suppressed 
and  the  plasma  temperature  decreases  toward  the  surface  of  the  liquid  helium.  The  cold  plasma  in 
a  steady  state  is  thus  produced. 

This  work  is  supported  by  Asian  Office  of  Aerospace  Research  and  Development. 
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Electrostatic  attractive  force  between  dust  particles  in  a  complex  plasma  with  ion  flow  is  studied.  It 
is  shown  that  there  is  an  attractive  force  between  a  pair  of  dust  particles  along  the  ion  flow  as  well 
as  perpendicular  to  the  ion  flow.  The  attractive  force  perpendicular  to  the  flow  results  from  a  release 
of  thermodynamic  free  energy  in  charged  fine  particles,  while  an  attractive  force  associated  with  the 
wake  potential  acts  on  a  pair  of  dust  particles  aligned  with  the  ion  flow.  Recent  experimental 
observation  of  the  sharp  boundary  of  a  void  in  a  complex  plasma  is  interpreted  as  a  result  of  the 
attractive  force.  ©  2005  American  Institute  of  Physics.  [DOI:  10.1063/1.1978467] 


A  void  structure,  characterized  by  a  dust-free  ordered 
state  with  a  sharp  boundary  in  a  dusty  plasma,  was  observed 
in  experiments  in  microgravity  conditions1  as  well  as  in 
ground-based  laboratory  plasmas.  A  theory  was  presented  to 
explain  the  void  formation  based  on  the  balance  of  the  elec¬ 
trostatic  force  and  the  ion  drag  force  acted  on  a  dust 
particle.3  Since  accumulated  observations  of  the  void  in  a 
complex  plasma  challenged  the  theory,  alternative  theories 
including  thermophoretic  force,2,4  phenomenological  theory,5 
and  nonlinear  theory6,7  have  been  presented.  Meanwhile,  the 
ion  drag  force  reevaluated  recently  by  taking  the  screened 
Coulomb  potential  into  account  is  revealed  to  be  much  stron- 
ger  than  that  estimated  earlier.  Most  recent  experimental 
observation  confirmed  the  role  of  the  ion  drag  force  for  the 
formation  of  the  void.9  Voids  have  been  observed  not  only  in 
a  complex  plasma,  but  also  in  colloidal  suspensions  where  an 
attractive  pair  interaction  was  considered.10,11  Like-charge 
attraction  is  a  recent  subject  in  a  plasma  as  well  as  in  a 
colloid.12,13  In  this  Letter,  we  consider  a  complex  plasma, 
where  macroparticles  are  embedded  in  a  fully  ionized  plasma 
and  the  system  is  thermodynamically  closed,  and  we  show 
that  a  like-charge  attraction  could  align  dust  particles  along 
the  equipotential  line  on  a  void  boundary  perpendicular  to 
the  ion  flow  in  the  complex  plasma,  while  the  ion  flow  will 
produce  the  wake  potential  behind  dust  particles  along  the 
direction  of  the  flow. 

We  consider  a  complex  plasma  in  which  N  dust  particles 
are  embedded  in  a  fully  ionized  collisionless  plasma.  We 
neglect  the  gravitational  effect,  as  is  the  case  in  microgravity 
experiments.  The  radius  of  a  spherical  dust  particle,  a ,  is 
assumed  to  be  much  smaller  than  the  Debye  shielding  length 
XD.  The  charge  density  of  the  complex  plasma  is  given  by 

N 

p(x)  =  Zde^j  8{x  -  Xj)  +  en^x)  -  ene(x ) ,  (1) 

7=1 

where  Zde  is  the  charge  of  a  dust  particle  located  at  x=x/-, 
nt(x)  is  ion  density,  and  ne(x)  is  electron  density.  A  charge 
neutrality  condition  is  given  by 


ne  ~  nt  =  NZd/V ,  (2) 

where  V  is  the  volume  occupied  by  a  complex  plasma,  nt 
=  fd3xni(x)/V  and  ne-  fd3xne(x)/V.  The  electrostatic  poten¬ 
tial  produced  by  N  dust  particles  at  a  position  x  at  time  t  is 
given  by 


2  <?'k  x  f  —  e-i(0t<f>(k,co), 

k  J  Br  277" 


(3) 


where  Br  indicates  the  Bromwich  integration  path  and 


7>(k,  m) 


i47TZ^  Y  -ik-Tj 
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(4) 


The  plasma  in  consideration  is  assumed  to  be  characterized 
by  a  dielectric  constant 


e(k,<u)  =  1  + 


* De 


J2J 


k2  («  -  k  •  Vy)2’ 


(5) 


where  \f  is  the  velocity  of  the  flowing  ions,  k2)e 
=4rrnee2 1  kT(„  a)2i=4Trnie2/mi  with  the  electron  temperature 
Te,  the  ion  mass  mh  and  the  Boltzmann  constant  k.  The  di¬ 
electric  constant  (5)  characterizes  the  plasma  where  the  ion 
flow  is  dominant  and  the  collisions  are  negligible.  We  note 
that  the  collective  modes  expressed  by  Eq.  (5)  damp  for  the 
wave  numbers  k>kDe.  The  poles  of  Eq.  (3)  with  cf)( k,co)  by 
Eq.  (4)  are  located  at  co=0  and  co=<u(k),  where  e[k,co(k)] 
=  0  and  to(k)  =  k  •  \f±  cok.  Here  we  introduced  cok 
=  (opik/  ^k2+kpe.  We  obtain  the  electrostatic  potential  as 


0(x)  =  0 o(x)  +  01  (X)  +  0w(X)> 


(6) 
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and  </>i(x)  is  the  potential  contribution  from  the  pole  o) 
=  a>(k)  which  characterizes  the  strong  interaction  of  a  dust 
particle  with  the  collective  plasma  mode.  The  potential 
c t>w(x )  is  known  to  be  responsible  for  producing  an  oscillat¬ 
ing  wake  potential  behind  each  dust  particle  along  the  ion 
flow.14-18  For  a  single  dust  particle  (N=  1),  the  contributions 
from  the  potential  </>0(x)  and  </>i(x)  were  shown  to  be  can¬ 
celed  by  the  part  of  the  potential  <fiw(x)  resulting  in  the  wake 
potential  existing  only  behind  the  isolated  dust  particle  in  the 
supersonic  ion  flow.14  The  wake  potential  is  the  source  of  the 
attractive  force  between  dust  particles  along  the  ion  flow. 
Because  of  such  an  attractive  force,  dust  particles  align  them¬ 
selves  along  the  direction  of  the  ion  flow.  We  now  show  that 
for  a  multiparticle  system  (N^2),  the  potential  <£0(x)  is  re¬ 
sponsible  to  produce  an  attractive  force  between  dust  par¬ 
ticles  perpendicular  to  the  ion  flow. 

The  electrostatic  energy  in  the  presence  of  a  supersonic 
ion  flow  is  given  by 
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v  <ft(x)|: 
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The  energy  depends  explicitly  on  Zd  and  the  coupling  param¬ 
eter  rj(=e2),  and  implicitly  on  ne  and  nt  through  Zd  and  kDe. 
We  note  that  the  electrostatic  energy  is  induced  by  a  pair 
displacement  x;--x7-,  a  novel  feature  of  a  multiparticle  sys¬ 
tem.  Such  a  displacement  necessarily  introduces  the  volume 
change  in  the  system  through  the  relationship  G=F+PV, 
where  G  is  the  Gibbs  free  energy  and  F  is  the  Helmholtz  free 
energy.  The  thermodynamic  excess  energy  (Gibbs  energy) 
thus  may  be  evaluated  as 


G=  2j  \  ns—  +  ns——  - —dr] 

s=e,i  '  ^ns  dns  vZdJ  Jq  7] 
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Here  we  use  ned(k2+k2D^~l  /  dne=-k2De(k2+k2D^)~2  and 
nedZd!  dne  +  riidZJ  rJEi  -  Zd.  The  k  summation  may  be  evalu¬ 
ated  through  the  integration 

y  f00  f 77  r2jr 

E  - >  7 — 77  k2dk  sin  0  dO  dtp, 
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where  k  is  chosen  to  be  along  the  z  axis,  6  and  cp  are  polar 
angles  of  x—Xy  with  respect  to  the  z  axis.  The  excess  energy 
becomes 


kDAj 


FIG.  1.  Normalized  potential  as  a  function  of  the  normalized  interparticle 
distance. 


c=Et/,, 

Uj 


Uij  = 


2Z^ 

Ru 


(10) 


where  Rij=\xi-xj\  is  the  interparticle  distance,  and  the  force 
acting  on  a  dust  pair  is  given  by 


dUjj 

F"~~  A, 


jzy 

rI 


1  +kDeRij--k2DeR2)e-k»A.  (11) 


Figure  1  shows  the  normalized  potential  as  a  function  of  the 
normalized  interparticle  distance.  The  potential  becomes 
minimum  at  Rtj = 4 . 83^  =  R0 . 

In  a  complex  plasma  experiment  with  parallel  electrode 
plates,  a  potential  has  its  maximum  at  the  center  of  the 
plasma  between  electrodes.  The  equipotential  lines  close 
themselves  making  circles  or  deformed  ellipsoidal  lines  in 
the  cross  section.19  Thus  the  electric  field  E  is  directing  out¬ 
ward  from  the  center  driving  the  electrons  inward,  while  pro¬ 
ducing  the  flux  of  ions  outward.  Negatively  charged  dust 
particles  are  pushed  inward  to  the  center  by  the  electric 
force,  FE=ZdeE,  while  the  ion  drag  force,  Fion 
=miJ\[(Tc(v)  +  (Ts(v)]vfi(\)d3v,  pushes  dust  particles  from 
the  center  toward  the  outer  region.9  Here  /)(v)  is  the  ion 
distribution  function,  and  the  collection  and  the  scattering 
cross  sections  are  given  by  ac(v)  =  7ra2(l-\-2b0/a),  crs 
=  47tZ?o  ln[(b0+\D)/(b0  +  a)~\,  respectively,  where  b0 
=Zde2/miv2.20  Due  to  the  balance  of  the  two  forces,  dust 
particles  characterized  by  the  charge  Zde  find  their  positions 
on  an  equipotential  line.  Some  deviation  of  positions  may  be 
expected  by  the  presence  of  other  forces  like  neutral  drag 
force,  gravitational  force,  and  Coulomb  force.21  Each  dust 
particle  on  the  equipotential  line  will  be  accompanied  by 
aligned  dust  particles  along  the  ion  flow  due  to  the  wake 
potential,  resulting  in  the  void  structure  in  the  central  part  of 
the  plasma.  Dust  particles  of  equal  charges  are  positioned  on 
the  same  equipotential  line  perpendicular  to  the  direction  of 
the  ion  flow  and  are  subject  to  the  force  given  by  Eq.  (11). 
Figure  2  describes  the  contour  plot  of  the  potential  around  a 
dust  particle  superimposed  on  the  wake  potential  along  the 
ion  flow  with  Mach  number  M= 2.  The  solid-black  three 
circles  are  dust  particles  on  an  equipotential  line,  while  the 
white  circular  band  around  a  dust  particle  in  the  middle 


Downloaded  24  Jul  2005  to  133.34.150.160.  Redistribution  subject  to  AIP  license  or  copyright,  see  http://pop.aip.org/pop/copyright.jsp 


070705-3  Attractive  force  on  like  charges  in  a  complex  plasma 


V 

EQUI>OTENTIAL 

LINE 


FIG.  2.  Contour  plot  of  the  potential  around  dust  particles  in  the  ion  flow. 
Three  dust  particles  are  shown  by  solid  black  circles.  Wake  potentials  are 
shown  in  the  direction  of  the  ion  flow  behind  the  particles  and  the  equipo- 
tential  lines  due  to  the  particle  in  the  middle  are  also  shown.  The  two 
neighboring  particles  are  located  on  the  potential  dips  created  by  the  middle 
particle. 

shows  the  presence  of  a  potential  minimum  with  a  radius  R0. 
The  neighboring  two  dust  particles  on  an  equipotential  line 
are  settled  at  the  minimum  potential  created  by  the  dust  par¬ 
ticle  in  the  center  with  the  interparticle  distance  of  R0.  The 
interparticle  distance  perpendicular  to  the  flow  is  indepen¬ 
dent  of  the  Mach  number,  while  the  supersonic  ion  flow  pro¬ 
duces  the  pattern  of  wake  potential.  With  a  smaller  Mach 
number,  the  distance  between  potential  minimums  along  the 
flow  behind  a  dust  particle  becomes  smaller.14  It  is  straight¬ 
forward  to  calculate  the  interparticle  distance  by  including 
the  finite  radius  a  of  a  dust  particle  in  a  dust  density  as 
Zde^Z  S(\x-xi\-a)/47ra2.  The  interparticle  distance  is  then 
given _ by  Rtj = ( 1  +  a  coth  kDea )  ( 1 

+  Jl+2/(l+kDea  coth  kDea)),  as  was  first  discussed  in  the 

22 

context  of  colloids. 

In  conclusion,  the  electrostatic  attractive  force  between 
like  charges  in  a  complex  plasma,  where  charged  macropar¬ 
ticles  are  embedded  in  a  fully  ionized  plasma,  is  shown  to  be 
responsible  for  the  sharp  boundary  formation  associated  with 
a  void  which  is  prompted  by  the  ion  drag  force.  The  attrac¬ 
tive  force  may  be  interpreted  as  a  release  of  thermodynamic 
free  energy  among  charged  fine  particles  in  the  presence  of 
ion  flow,  while  the  attractive  force  associated  with  the  wake 
potential  may  be  interpreted  as  a  result  of  phonon  exchange 
in  the  presence  of  ion  flow  in  a  complex  plasma.16  The  mea¬ 
surement  of  interparticle  distances  may  be  used  to  estimate 
plasma  parameters.  The  interparticle  distance  perpendicular 
to  the  boundary  will  provide  the  ion  flow  velocity,  while  the 
interparticle  distance  along  the  boundary  will  provide  the 
Debye  screening  length.  Although  the  exact  plasma  param¬ 
eters  at  the  boundary  of  the  void  may  be  difficult  to  measure, 
some  typical  interparticle  distance  of  small  grains  with  ra¬ 
dius  on  the  order  of  a  few  micrometers  has  been  reported  to 
be  on  the  order  of  0.1  mm  in  the  range  of  a  few  ion  Debye 
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lengths  and  less  than  an  electron  Debye  length  in  micrograv¬ 
ity  conditions  with  ne^ni  —  109  cm-3,  Te~  3  eV,  and  Tt 
~  0.03  eV.  It  should  be  noted  that  the  proposed  void  struc¬ 
ture  is  operative  in  a  complex  plasma  which  is  thermody¬ 
namically  closed.  In  our  analysis  we  have  neglected,  as  is 
shown  by  the  dielectric  constant  (5),  the  effect  of  ion-neutral 
collisions  and  thermal  spread  in  ion  velocities  which  are 
commonly  present  in  laboratory  experiments.  Since  a  weakly 
ionized  plasma  with  neutrals  and  finite  ion  temperature  is 
characterized  by  an  open  dissipative  system  in  dust  plasma 
experiments,  some  deviation  on  the  interparticle  distance  for 
attractive  force  to  occur  may  be  expected. 

This  research  was  supported  by  Japan  Society  for  the 
Promotion  of  Science  Grant-in-Aid  for  Scientific  Research 
(B)  under  Grant  No.  (2)  1  6  3  4  0  1  7  9. 

!G.  E.  Morfill,  H.  M.  Thomas,  U.  Konopka,  H.  Rothermel,  M.  Zuzic,  A. 
Ivlev,  and  J.  Goree,  Phys.  Rev.  Lett.  83,  1598  (1999);  A.  P.  Nefedov,  G. 

E.  Morfill,  V.  E.  Fortov,  H.  M.  Thomas,  H.  Rothermel,  T.  Hagl,  A.  V. 
Ivlev,  M.  Zuzic,  B.  A.  Klumov,  A.  M.  Lipaev,  V.  I.  Molotkov,  O.  F. 
Petrov,  Y.  P.  Gidzenko,  S.  K.  Krikalev,  W.  Shepherd,  A.  I.  Ivanov,  M. 
Roth,  H.  Binnenbruck,  J.  A.  Goree,  and  Y.  P.  Semenov,  New  J.  Phys.  5, 
33.1  (2003);  Y.  Hayashi,  Jpn.  J.  Appl.  Phys.,  Part  1  44,  1436  (2005). 

2G.  Praburam  and  J.  Goree,  Phys.  Plasmas  3,  1212  (1996);  D.  Samsonov 
and  J.  Goree,  Phys.  Rev.  E  59,  1047  (1999);  H.  Rothermel,  T.  Hagl,  G.  E. 
Morfill,  M.  H.  Thoma,  and  H.  M.  Thomas,  Phys.  Rev.  Lett.  89,  175001 
(2002). 

3R.  P.  Dahiya,  G.  V.  Paeva,  W.  W.  Stoffels,  E.  Stoffels,  G.  M.  W.  Kroesen, 
K.  Avinash,  and  A.  Bhattacharjee,  Phys.  Rev.  Lett.  89,  125001  (2002). 
4J.  Goree,  G.  E.  Morfill,  V.  N.  Tsytovich,  and  S.  V.  Vladimirov,  Phys.  Rev. 
E  59,  7055  (1999). 

5K.  Avinash,  Phys.  Plasmas  8,  2601  (2001). 

6K.  Avinash,  A.  Bhattacharjee,  and  S.  Hu,  Phys.  Rev.  Lett.  90,  075001 
(2003). 

7D.  Jovanovic  and  P.  K.  Shukla,  Phys.  Lett.  A  308,  369  (2003). 

8S.  A.  Khrapak,  A.  V.  Ivlev,  G.  E.  Morfill,  and  S.  K.  Zhdanov,  Phys.  Rev. 
Lett.  90,  225002  (2003). 

9M.  Kretschmer,  S.  A.  Khrapak,  S.  K.  Zhdanov,  H.  M.  Thomas,  G.  E. 
Morfill,  V.  E.  Fortov,  A.  M.  Lipaev,  V.  I.  Molotkov,  A.  Ivanov,  and  M.  V. 
Turin,  Phys.  Rev.  E  71,  056401  (2005). 

10D.  G.  Grier  and  J.  C.  Crocker,  Phys.  Rev.  E  61,  980  (2000). 
nN.  Ise,  Proc.  Jpn.  Acad.,  Ser.  B:  Phys.  Biol.  Sci.  78,  129  (2002). 

12V.  N.  Tsytovich,  Comments  Plasma  Phys.  Controlled  Fusion  15,  349 
(1994);  V.  N.  Tsytovich  and  G.  E.  Morfill,  Plasma  Phys.  Rep.  28,  171 
(2002). 

13T.  M.  Squires  and  M.  P.  Brenner,  Phys.  Rev.  Lett.  85,  4976  (2000). 

14 O.  Ishihara  and  S.  V.  Vladimirov,  Phys.  Plasmas  4,  69  (1997). 

15M.  Nambu,  S.  V.  Vladimirov,  and  P.  K.  Shukla,  Phys.  Lett.  A  203,  40 
(1995);  S.  V.  Vladimirov  and  O.  Ishihara,  Phys.  Plasmas  3,  444  (1996). 
160.  Ishihara  and  S.  V.  Vladimirov,  Phys.  Rev.  E  57,  3392  (1998). 

17D.  Winske,  W.  Daughton,  D.  S.  Lemons,  and  M.  S.  Murillo,  Phys.  Plasmas 
7,  2320  (2000);  M.  Lampe,  G.  Joyce,  G.  Ganguli,  and  V.  Gavrishchaka, 
ibid.  7,  3851  (2000). 

18K.  Takahashi,  T.  Oishi,  K.  Shimonai,  Y.  Hayashi,  and  S.  Nishino,  Phys. 
Rev.  E  58,  7805  (1998). 

19M.  R.  Akdim  and  W.  J.  Goedheer,  Phys.  Rev.  E  65,  015401(R)  (2001). 
20S.  A.  Khrapak,  A.  V.  Ivlev,  G.  E.  Morfill,  and  H.  M.  Thomas,  Phys.  Rev. 
E  66,  046414  (2002). 

210.  Ishihara,  Phys.  Plasmas  5,  357  (1998). 

22I.  Sogami  and  N.  Ise,  J.  Chem.  Phys.  81,  6320  (1984). 

23 G.  E.  Morfill,  M.  Rubin-Zuzic,  H.  Rothermel,  A.  V.  Ivlev,  B.  A.  Klumov, 
H.  M.  Thomas,  U.  Konopka,  and  V.  Steinberg,  Phys.  Rev.  Lett.  92, 
175004  (2004);  V.  V.  Yaroshenko,  B.  M.  Annaratone,  S.  A.  Khrapak,  H. 
M.  Thomas,  G.  E.  Morfill,  V.  E.  Fortov,  A.  M.  Lipaev,  V.  I.  Molotkov,  O. 

F.  Petrov,  A.  I.  Ivanov,  and  M.  V.  Turin,  Phys.  Rev.  E  69,  066401  (2004). 


Downloaded  24  Jul  2005  to  133.34.150.160.  Redistribution  subject  to  AIP  license  or  copyright,  see  http://pop.aip.org/pop/copyright.jsp 


